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Abstract: 

B — > ^ev transitions have recently been studied in the framework of QCD factorization. 
The attractiveness of this channel for such an analysis lies in the fact that, at least in 
the heavy quark limit, the only hadron involved is the B meson itself, so one expects 
a very simple description of the form factor in terms of a convolution of the B meson 
distribution amplitude with a perturbative kernel. This description, however, does not 
include contributions suppressed by powers of the b quark mass. In this letter, we cal- 
culate corrections to the factorized expression which are induced by the "soft" hadronic 
component of the photon. We demonstrate that the power-suppression of these terms is 
numerically not effective for physical values of the b quark mass and that they increase the 
form factor by about 30% at zero momentum transfer. We also derive a sum rule for As, 
the first negative moment of the B meson distribution amplitude, and find Xb = 0.6 GeV 
(to leading order in QCD). 
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1 Introduction 



With the B factories BaBar and Belle running "full steam", B physics has entered the 
era of precision measurements. The quality and precision of experimental data calls for 
a corresponding match in theoretical precision, in particular with regard to the analysis 
of nonleptonic decays. Notwithstanding the fact that a complete solution of the problem 
appears as elusive as ever, progress has been made in the description of nonleptonic B 
decays, in the heavy quark limit, which were shown to be amenable to perturbative QCD 
(pQCD) factorization [1, 2]. In this framework, B decay amplitudes are decomposed into 
a "soft" part, for instance a weak decay form factor, or, in general, an intrinsically non- 
perturbative quantity that evades further breakdown into factorizable components, and a 
"hard" part which can be neatly described in a factorized form in terms of a convolution 
of a hard perturbative kernel, depending only on collinear momenta, with one or more 
hadron distribution amplitudes that describe the (collinear) momentum distribution of 
partons inside the hadron. Factorization was shown to hold, for certain B decay channels, 
to all orders in perturbation theory [3], but to date could not be extended to include con- 
tributions that are suppressed by powers of the b quark mass. For the important channel 
B — > Kn, in particular, it was found that factorization breaks down for one specific class 
of power-suppressed corrections which are expected to be numerically relevant [2]. Reli- 
able alternative methods for calculating QCD effects in nonleptonic B decays are scarce, 
and only little is known about the generic size of power-suppressed corrections for instance 
in B — > nil [4], one of the "benchmark" channels for measuring the angle a of the CKM 
unitarity triangle. There is, however, one channel that can be treated both in pQCD 
factorization and by an alternative method, QCD sum rules: the B — > jivg transition. It 
has been shown recently [5, 6] that B — > 7 is indeed accessible to collinear factorization, 
in contrast to previous findings [7] which indicated the need to include also transverse 
degrees of freedom in the convolution. On the other hand, the I? — > 7 transition can also 
been investigated in the framework of QCD sum rules on the light-cone [8, 9]. The crucial 
point here is that the photon is not treated as an exactly pointlike object with standard 
EM couplings, but that, in addition to that "hard" component, it also features a "soft" 
hadronic component whose contribution to the decay amplitude must not be neglected. 
The "soft" component is related to the probability of a real photon to dissociate, at small 
transverse separation, into partons and resembles in many ways a (massless) transversely 
polarized vector meson; like it, it can be described by a Fock-state expansion in terms of 
distribution amplitudes of increasing twist. An analysis of these distribution amplitudes, 
including terms up to twist-4, has recently been completed [10] and comes in handy for 
an update and extension of the previous QCD sum rules analyses of B — > 7. Such a 
reanalysis is the subject of this letter and we include in particular one-loop radiative cor- 
rections to the contribution of leading twist-2 distribution amplitude to B — > 7. In the 
framework of pQCD factorization, the "soft" component of the photon leads to formally 
power-suppressed contributions and is hence neglected in Refs. [5, 6, 7]. As we shall show 
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in this letter, this suppression is not effective numerically due to the large value of the 
matrix-element governing its strength: the magnetic susceptibility of the quark conden- 
sate. Our findings indicate that such contributions are likely to be nonneglible also in 
other channels involving photon emission, notably B — > K*^ and B — > p7, which are 
treated in [11, 12]. We also compare the QCD sum rule for the hard part of the 
amplitude with the pQCD result and derive a sum rule for \ B , the first negative moment 
of the B meson distribution amplitude. 

2 Definition of Relevant Quantities and Outline of 
Calculation 

First of all, let us define the form factors that describe the B — > 7 transition: 

-^=(7(e*, ? )|fZ7 M (l-7 5 )fe|S-(p B )) = -F v e flupa e*"v^ + tF A [e;(vq)-q II (e*-v)}, (1) 

where e* and q are the polarisation and momentum vector of the photon, respectively, 
and v = ps/fTiB is the four-velocity of the B meson. The definition of the form factors 
Fa,v in (1) is exactly the same as in Ref. [5]. 1 The form factors depend on p 2 = (p B — q) 2 
or, equivalently, on the photon energy E 1 = (m 2 B — p 2 )/(2m B ). p 2 varies in the physical 
region < p 2 < m 2 B , which corresponds to < E 1 < m B /2. 

The starting point for the calculation of the form factors from QCD sum rules is the 
correlation function 

n>,g) = i I rf 4 a:e^^^(7(e*,g)|T{«(a:)7^1-75)K^(0)^(0)}|0) 

= -IV W e* W + iK A [el{p ■ q) - g M (e* -p)] + ... (2) 

The dots refer to contact terms which appear for pointlike photons and for a discussion 
of which we refer to Ref. [13]; the treatment of the soft component of the photon involves 
nonlocal operators, as we shall discuss below, and gauge-invariance of I1 M is realized 
explicitly, without contact terms, by working in the background field method. 

The method of QCD sum rules [14] exploits the fact that the correlation function 
contains information on the form factors in question: expressing Ily(A) via a dispersion 
relation, one has 

f B m B F v{A) r°° ds 2 
Uv ^ = m b (m 2 B - P 2 B ) + L S—f B PV ^ P ^ (3) 

where the first term on the r.h.s. is the contribution of the ground state B meson to the 
correlation function, featuring the form factors we want to calculate, and the second term 

1 Thc difference in sign in the 1st term on the r.h.s. is due to our conventions for the epsilon tensor: 
we define Tr [7^7,^ 7^75] = Aie^ upa , in contrast to [5], where a different sign was chosen. 
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Figure 1: Diagrams with hard photon emission. 

includes all other states coupling to the pseudoscalar current bi'y^u, above the threshold 
so- /b is the decay constant of the B and defined as {B\bi^u\fy = /Bmg/mj. The form 
factors F A y are obtained, in principle, by equating the dispersion relation (3) to Ha,v 
calculated for Euclidean p 2 B by means of an operator product expansion. The sum over 
higher states, the 2nd term on the r.h.s. of (3), is evaluated using quark-hadron duality, 
which means that the hadronic spectral density pv(A) is replaced by its perturbative 
equivalent. In order to reduce the model-dependence associated with that procedure, one 
subjects the whole expression to a Borel-transformation, which results in an exponential 
suppression of the continuum of higher states: 

B M ,n v{A) = f* m * F v^ e - ml/ M> + r dspV « {s)e -s /M ^ (4) 

nib Js 

The relevant parameters of the sum rule are then M 2 , the Borel parameter, and s , the 
continuum threshold, and our results for the form factors will depend (moderately) on 
these parameters. 

As already mentioned in the introduction, the photon does not only have EM pointlike 
couplings to the quarks, but also soft nonlocal ones, and we write F A ,v — F A a ^ + ^av- 
The separation between these two components is of course scheme-dependent and we 
define them in the MS scheme. For the pointlike component of the photon, IT M is just 
a three-point correlation function, described by the triangle diagrams shown in Fig. 1. 
These diagrams, and also the leading nonperturbative correction which is proportional to 
the quark condensate, have been calculated in Refs. [8, 9], and we confirm the results; 
we will come back to the hard contributions in Sec. 4. The calculation of radiative 
corrections to the three-point function, although highly desirable, is beyond the scope 
of this letter whose main emphasis is on calculating the soft photon contributions which 
become relevant for a certain configuration of virtualities, namely ml — p%> 0(Aqcd?™&) 
and ml — p 2 > 0(AQ CD m b ), i.e. ^> Aq CD . In this regime, the integral in Eq. (3) is 
dominated by light-like distances and can be expanded around the light-cone: 

U v{A) (p 2 B ,E 1 )=Q uX M(uu)(fi F )J2 f du^ n \u ] p F )T^ A) (u ] p 2 B ,E l] p F ). (5) 

This is a factorization formula expressing the correlation function as a convolution of 



3 



genuinely nonperturbative and universal distribution amplitudes (DA) with process- 
dependent hard kernels T^ n \ to be calculated in perturbation theory; the overall factor 
QuX{uu) (Q u = 2/3) is pulled out for later convenience. In the above equation, n labels 
the twist of operators and n F is the factorization scale. The restriction on E 1 implies 
that F^y cannot be calculated for all photon energies; to be specific, we restrict ourselves 
to E 1 > 1 GeV. The factorization formula holds if the infrared divergencies occurring in 
the calculation of are such that they can be absorbed into the universal distribution 
amplitudes (j)^ and if the convolution integral converges. We find that this is indeed the 
case, at least for the leading twist- 2 contribution and to first order in a s . We also would 
like to stress that the above factorization formula has got nothing to do with the pQCD 
factorized expression for F^ v d derived in [5, 6] (although we will discuss the relevance of 
our findings as compared to theirs in Sec. 4) and that it is valid for arbitrary values of 
the b quark mass — there is no need to restrict oneself to the heavy quark limit. 

Let us now define the photon DAs that enter Eq. (5). As mentioned above, we shall 
work in the background field gauge, which means that the expression (3), formulated for 
an outgoing photon, gets replaced by a correlation function of the time-ordered product 
of two currents in the vacuum, which is populated by an arbitrary EM field configuration 

Fa/3, 

n F (p,q) = i I d A xe^^={Q\T{u{x)r{l - lb )b(x)b(ti)i lb u{<S)}\()) F , (6) 
J y Ana 

where the subscript F indicates that an EM background field is included in the action. 
The calculation is explicitly gauge-invariant, and it is only in the final step that we select 
one specific field configuration corresponding to an outgoing photon with momentum q: 

F a fs(x) -> -i(e* a q p - e* p q a )e iqx . 

Following Ref. [10], we define the leading-twist DA as the vacuum expectation value of 
the nonlocal quark-antiquark operator with light-like separations, in the EM background 
field configuration F a p: 

l 

(0\q{z)<J a p[z,-z] F q(-z)\0) F = e q x{qq) j "du F a/3 (-^)^ 7 (u) , (7) 

o 

where z 2 = and [z, —z] F is the path-ordered gauge-factor 

[z, -z] F = Pexp ^dt2z^gA^(2t - l)z) + e q B<*((2t - l)z)]} , (8) 

including both the gluon field = A^\ a /2 and the EM background field B^. e q is 
the electric charge of the light quark q, e.g. e u = 2/3^/Aira, (qq) the quark condensate 
and x its magnetic susceptibility - defined via the local matrix element (0\qa al 3q\0) F = 
e qX (QQ.) F a /3, which implies the normalization Jq du(f>^{u) = 1. 7 can be expanded in 
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Figure 2: Leading twist-2 soft photon contributions to the correlation function (2). 



terms of contributions of increasing conformal spin (cf. [15] for a detailed discussion of 
the conformal expansion of DAs) : 



(u, //) = 6u(l 



u 



i+ J2 Mv)c 3 n /2 (2u-i) 

n=2,4,... 



(9) 



where C^J 2 are Gegenbauer polynomials and < u < 1 is the collinear momentum fraction 
carried by the quark. The usefulness of the conformal expansion lies in the fact that the 
Gegenbauer moments n (/-O renormalize multiplicatively in LO perturbation theory. At 
higher twist, there exists a full plethora of photon DAs, which we refrain from defining in 
full detail, but refer the reader to the discussion in Ref. [10], in particular Sec. 4.1. 

The task is now to calculate the hard kernels T^ n \ to 0(a s ) in leading-twist and 
tree-level for higher twist. As for the former, the relevant diagrams are shown in Fig. 2. 
The diagrams are IR divergent, which divergence can be absorbed into the DA x(mi)(/> 7 ; 
we have checked that this is indeed the case. A critical point in pQCD calculations 
involving heavy particles is the possibility of soft divergencies, which manifest themselves 
as divergence of the u integration at the endpoints; the success of the approach advocated 
in [1, 2] relies precisely on the fact that it yields a factorization formula for nonleptonic 
decays where these divergencies are absent in the heavy quark limit (but come back at 
order 1/mb). We find that there are no soft divergencies in our case 2 , which reiterates 
what has been found for other heavy-light transitions, cf. [16, 17, 18]. We thus confirm 
the factorization formula (5) to 0(a s ) at twist-2. Explicit expressions for the spectral 

2 Note that this statement applies to the full correlation function with finite nib, not only to the heavy 
quark limit. 



densities of the diagrams are too bulky to be given here; they can be obtained from the 
authors. At this point we only note that, to leading-twist accuracy, II a = U v . 
As for the higher-twist contributions, we obtain the following results: 



f 1 ip^iu) 1 _ f 

ILi = -Qufz-t / du — m b + -Q u (uu) 

Jo s z 2 Jo 



1 + 



2m 2 



(10) 



n 



V 



-Q u {uu) f dv f Va^-(1 - 2v) + \q u {uu) f dv f Va^- 
-2Q u (uu) f dv f Va^- [Ti(a) - (1 - 2v)T 2 (a) + (1 - 2v)T 3 (, 

1 ,1 

T^Quhi / du — m h 

Zi jo s 

-Qu{uu) f dv f Va^f- + \q u {uu) f dv f Va^-{1 - 2v) 

J J S J J s 



(11) 



p-q 



-2Q u (uu) / dv / Va^— i 71(a) - T 2 (a) + (1 - 2t;)T 3 (a) - (1 - 2v)T 4 (a) 



to J s° 

where s = ml — {p + -ug) 2 and s = m 2 — (p + (q) 2 with £ = (a^ — a q + {2v — l)a g + l)/2. 
Definitions and explicit expressions for the higher twist DAs ^ v \u) and ^ a \u) (two- 
particle twist-3), h(u) (two-particle twist-4), S(a), S(a) (three-particle twist-3) and Ti(a) 
(three-particle twist-4) can be found in [10]. Note that in contrast to the leading-twist 

results, n^ ighertwis V n^ ighcrtwist . 



3 Numerical Results for F^ v 

Before presenting numerical values for the soft part of the form factors, we first discuss 
the numerical input to the sum rule (4). As for the photon DA, we use the values 
and parametrizations derived in [10], notably for the normalization of the twist-2 matrix 
element (7): 

(x(uu))(l GeV) = -(0.050 ± 0.015) GeV. 

As discussed in [10], there is no conclusive evidence for 7 to differ significantly from its 
asymptotic form, so we set 

7 (w) = 6u(l — u). 

The remaining hadronic matrix elements characterizing higher-twist DAs are detailed in 
[10]. Note that we evaluate scale-dependent quantities at the factorization scale [i 2 F = 
m 2 B — m 2 [19]; the dependence of the form factors on fip is very small, as all numerically 
sizeable contributions are now available to NLO in QCD, which ensures good cancellation 
of the residual scale dependence. 
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As for the remaining parameters occurring in (4), we have the sum rule specific param- 
eters M 2 and s , that is the Borel parameter and the continuum threshold, respectively. 
In addition, the sum rule depends on mj, the b quark mass, and the leptonic decay 
constant of the B. /# can in principle be measured from the decay B — > £vt, which, due to 
the expected smallness of its branching ratio, BR~ O(10 -6 ), has, up to now, escaped ex- 
perimental detection. is one of the best-studied observables in lattice-simulations with 
heavy quarks; the current world-average from unquenched calculations with two dynami- 
cal quarks is fs = (200 ±30) MeV [20]. It can also be calculated from QCD sum rules: the 
most recent determinations [21] include 0(a 2 ) corrections and find (206 ± 20) MeV and 
(197 ±23) MeV, respectively. For consistency, we do not use these results, but replace f B 
in (4) by its QCD sum rule to 0(a s ) accuracy, including the dependence on s and M 2 , 
and use the corresponding "optimum" ranges of continuum threshold and Borel param- 
eter also in evaluating the Borel-transformed correlation function T1v(a), i-e. the l.h.s. of 
(4). For the b quark mass, we use an average over recent determinations of the MS mass, 
m b Ms(jn b ) = (4.22 ± 0.08) GeV [22, 23], which corresponds to the one-loop pole-mass 
^iL-poic = (4.60 ± 0.09) GeV. With these values we find f B = (192 ± 22) GeV (the error 
only includes variation with m b and M 2 , at optimized sq), in very good agreement with 
both lattice and QCD sum rules to 0(a 2 ) accuracy. For m b = (4.51,4.60,4.69) GeV the 
optimized s are (34.5, 34.0, 33.5) GeV 2 , and the relevant range of M 2 is (4.5-8) GeV 2 . 

In Fig. 3 we plot the different contributions to Fy°^(0) as function of M 2 , for m b = 
4.6 GeV, So — 34 GeV 2 and the central value of x(uu). It is evident that the sum rule is 
dominated by twist- 2 contributions and that both radiative corrections and higher-twist 
terms are well under control. Note also the minimal sensitivity to M 2 which indicates 
a "well-behaved" sum rule. Varying m^, s and the other input parameters within the 
ranges specified above, we find 

i^ oft (0) = 0.07 ± 0.02, Ff oft (0) = 0.09 ± 0.02. (12) 

As mentioned before, the above results are obtained using the asymptotic form of the 
twist-2 photon DA. Although there is presently no evidence for nonzero values of higher 
Gegenbauer moments, the (fi n in (9), it may be illustrative to estimate their possible 
impact on the form factors. As a guideline for numerics, we choose 4>1 to be equal to (p^, 
its analogue for the transversely polarized p meson, as determined in [24]: 02 X (1 GeV) = 
0.2 ± 0.1. In Fig. 4 we plot the twist-2 contribution to the form factors obtained with the 
asymptotic 7 and the corrections induced by nonzero ^{jXp)- It is clear that the effect 
is moderate and at most about 20%. It is also interesting to note that positive values of 
02; which are in accordance with assuming p meson dominance for the photon, increase 
the form factor, which means that our results are likely to be an underestimate of the soft 
contributions rather than the contrary. 

In Fig. 5, we show the dependence of F^y(p 2 ) on the momentum transfer p 2 , including 
the variation of all input parameters in their respective ranges. The form factors can be 
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Figure 3: The soft photon contributions to Fy (a)(0) in dependence of the Borel parameter M 2 . 
The twist-2 contributions to both form factors are equal. The plot refers to mi = 4.6 GeV and 
s = 34 GeV 2 . 
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Figure 4: Contributions to Fy^ from different values (f>2 of the 2nd Gegenbauer moment of 
the photon DA (m b = 4.6 GeV, s = 34 GeV 2 ). 
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Figure 5: Dependence of Fy^ on the momentum transfer p 2 . Solid line: light-cone sum rule 
for central values of input parameters. Dashed lines: spread of Fy^ upon variation of input 
parameters. The exact parameter sets for the curves (a) to (c) are listed in Tab. 1. Note that 
the uncertainty originating from M 2 and sq is less than 10%. 
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Table 1: Input parameter sets for Fig. 5 and fit parameters for Eq. (13). 



fitted by the following formula: 

nup 2 ) = - ; ,:z>:~: — ■ m 



1 - a A ( V ) (p 2 /m 2 B ) + b A (v) {p 2 /m 2 B) 



The fit parameters F A °^rJ0), cla(v) an d &A(v) for each curve in the figure are given in Tab. 1. 
The above formula fits the full sum rule results to 1% accuracy for < p 2 < 17GeV 2 , 
which corresponds to 1 GeV < E 1 < ms/^. Note that the uncertainty induced by M 2 and 
s is very small: up to 5% for F v and up to 10% for Fa- The main theoretical uncertainty 
comes from x{qq)- 



4 Calculation of A# - Comparison to pQCD 

As mentioned above, the B — > 7 form factors have also been calculated in the framework 
of pQCD factorization [7, 5, 6]. This approach employs the limit m h — > 00 (HQL), in 
which the soft contributions vanish. The form factors F A ^ HQL are equal and at tree 
level given by 

V 7 Q + 7 73 

where the light-cone DA <3>+ of the -B meson depends on the momentum of the light 
spectator quark, k + . The calculation of radiative corrections to this formula has been the 
subject of a certain controversy, cf. Refs. [7, 5, 6]. The parameter A73 does not scale with 
771^ in the HQL and hence is of natural size O(Aqcd); it has been quoted as A# = 0.3 GeV 
[1] and \b = (0.35 ± 0.15) GeV [5], but without calculation. Although presently any 
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statement about the numerical size of \ B appears slightly precarious, since Q B and hence 
\ B depend in a yet unknown way on the factorization scale /ip, we nonetheless venture 
to present the (to the best of our knowledge) first calculation of X B . To that purpose, we 
recall that the hard contribution to Fa,v can be obtained, in the QCD sum rule approach, 
from the local contributions to the correlation function IT^y, Eq. (2), which, to leading 
order in perturbation theory, correspond to the diagrams shown in Fig. 1 and have been 
calculated in [8]. In order to extract X B via (14) from the local QCD sum rule for F^ffi, 
we first have to find its HQL. QCD sum rules in the heavy quark limit have actually been 
studied in quite some detail, cf. [25], with the following result for the scaling relations of 
the sum rule specific parameters M 2 and s : 

M 2 — > 2rribT, so — > m 2 + 2m b uj§. (15) 

Applying these relations to the sum rule for Fy^ d and using (14), we obtain the following 
expression: 

e-^ f -^^ = 4-r d uue-^. (16) 
m fe £/ 7 \ B i\ l E 1 Jo 

We note that the factor 1/E 1 on the r.h.s. arises automatically in the HQL of the correla- 
tion function. A is the binding energy of the b quark in the B meson, A = m B — mj. On 
the l.h.s. of (16), the expression f B m 2 B /m h still contains l/m b corrections; in the rigorous 
HQL we have f B m 2 B /m b — > / 2 tat . / stat is known both from lattice calculations and QCD 
sum rules; in the same spirit that we applied for calculating F^y, we replace, for the 
numerical evaluation of As, f 2 tat by its sum rule [25]: 

fLe~^ = - 2 Tdu^e-^, (17) 



7T 2 JO 



where we have suppressed (small) condensate contributions. Combining (16) and (17), 
we find 

Ab = -7^5 — ■ (18) 

/ duue- u/T 
Jo 

This is our sum rule for X B , which, admittedly, is the first rather than the last word in the 
story of how to calculate X B . It can and should be improved by including both radiative 
(which will also settle the issue of scale dependence) and nonperturbative corrections. For 
fixing the sum rule parameters c^o and r, we exploit the fact that A = m B — m h pa 0.7 GeV 
is known; the corresponding sum rule reads 

I"* duu 3 e-" /T 

A = Jj h* r— > ( 19 ) 

duu 2 e-"l T 



Jo 



which can be derived from (17) by taking one derivative in 1/r. None of the sum rules 
(16), (17), (18) is "well-behaved" in the sense that they feature no stability plateau in r, 
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since the perturbative term is not counterbalanced by a nonperturbative one. Nonetheless, 
taking 0.5 GeV < r < 1 GeV as indicated by our preferred range of M 2 used above and 
the scaling laws (15), we find lo — 1 GeV and then, from (18), 

\ B ~ 0.57 GeV. 

Given the neglect of 0(a s ) corrections and nonperturbative terms, it is difficult to at- 
tribute an error to that number. For this reason, we check if our result is compatible with 
the local duality approximation, which corresponds to the limit r — > oo. From (16) and 
(17) we then obtain the interesting relation 




which for m b = 4.6 GeV implies A^ = 0.6 GeV. We thus conclude that the value of X B is 
set by A rather than Aqcd and depends strongly on the actual value of m&; at present, 
nothing meaningful can be said about the error associated with X B and we thus quote as 
our final result 

A B = 0.6GeV. (20) 

We are now in a position to compare the numerical size of the pQCD result F^ffA H q L to 
Fffiy At tree level, we have Fj^? ))HQL (£7 7 = m B /2) = 0.22, according to (14). Notwith- 
standing the additional hadronic uncertainties involved at NLO in QCD, we employ the 
models for $^ advocated in [5] to obtain ^^'hqlC-^y = m #/2) — 0.21 ± 0.01 for 
[/,% = m 2 B — mf, our choice of the factorization scale. This number has to be compared to 
(12), the soft contributions. We conclude that F^y/F^^ HQL « 0.3 at maximum photon 
energy so that the parametrical scaling F^y ~ 0(l/m b ) is numerically relaxed. 



5 Summary and Conclusions 

The relevance of B physics for extracting information on weak interaction parameters and 
new physics is limited by our lack of knowledge on nonperturbative QCD. Recent progress 
in describing the notoriously difficult nonleptonic decays in perturbative QCD factoriza- 
tion has raised hopes that a sufficiently accurate solution to the problem is around the 
corner. As much as this is a highly desirable goal, it is nonetheless necessary to critically 
examine the theoretical uncertainty of the method, which, at least at present, is set by 
the restriction to the heavy quark limit. A direct theoretical test of pQCD factorization 
in nonleptonic decays is currently not feasible, and any significant experimental deviation 
of measured decay rates or CP asymmetries from pQCD predictions is as likely to be 
attributed to new physics effects as to uncertainties in the predictions themselves. An 
indirect theoretical test becomes however possible in the admittedly phenomenologically 
not very attractive channel B — > '-fev, where alternative methods of calculation exist and 
allow one to assess effects suppressed by powers of m&. In this letter, we have calculated 
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corrections to the hard pQCD form factors, which are parametrically suppressed by one 
power of m b . These "soft" corrections are induced by photon emission at large distances 
and involve the hadronic structure of the photon. We have also presented the first cal- 
culation of Xb, the first negative moment of the B meson distribution amplitude, a very 
relevant parameter for pQCD calculations. The calculation is admittedly rather crude, 
but amenable to improvement. Comparing the numerical size of the pQCD result to 
the soft contributions, we found that the latter are indeed sizeable, 0(30%). This result 
implies an immediate caveat for pQCD analyses involving photon emission, in particular 
B — > K*y and B — > py, e.g. [11, 12]. In a wider sense, it also adds a possible question- 
mark to pQCD analyses of purely hadronic B decays and emphasizes the relevance of 
power-suppressed corrections to the heavy quark limit. 
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